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Abstract 

The harmonic inversion method is apphed in the case of the hydrogen atom 
in a magnetic field to extract classical information from the quantum photo- 
ionization cross-section. The study is made close to a saddle-node bifurcation 
for which the usual semi-classical formulas give diverging contributions. All 
quantities (actions, stabilities and Maslov indices) for real orbits above the 
bifurcation and for complex ghost orbits below the bifurcation, are found to 
be in excellent agreement with the modified semi-classical predictions based 
on a normal form approach. 

PACS number(s): 31.15.Gy, 05.45.Mt, 03.65.Sq, 32.60.i 
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I. INTRODUCTION 



The hydrogen atom in a strong magnetic field is one of the most appeahng system to 
study chaotic effects in quantum mechanics: it has the minimum number of degrees of 
freedom, all quantum and classical properties can be computed in a exact way [0,0. Some 
other systems share these properties, but the difference is that the chaotic regime for the 
magnetized hydrogen atom is also observed experimentally which is the case of only few 
systems, as helium atom , hydrogen atom in a microwave field , resonant tunneling 



diode pHllI, experiments with cold atoms |12| , p!3| . Chaos being defined at the classical 



level as exponential sensitivity on initial conditions, a major step in the understanding of 
chaos at the quantum level has been made trough the Gutzwiller trace formula, linking 
the quantum density of states to the classical periodic orbits of the system [Q. Later 
on, similar formula were derived to explain experimental observation, like photo-ionization 
cross-section, in terms of classical orbits (closed orbits in this case) |T3,TB]. 

The agreement between the numerical and experimental quantum results was found to be 
good, except for parameters (energy, magnetic field) too close to bifurcation points, at which 
the semi-classical amplitudes are diverging. However, these divergences have been later 
understood and classified using a normal form approach fl,p!7|. Especially, it emphasized 



the contribution of complex "ghost" |]TH[ orbits (where both position and momentum are 
made complex quantities) in the photo-ionization cross- sect ion. This again was found to 
be in a qualitatively good agreement with numerical and experimental data. Still, the 
comparison between the quantum results and the classical predictions being obtained by 
Fourier transforms of the quantum data, the accuracy is limited by the length of the available 
spectra (experimental or theoretical). For example, the imaginary part of the action of the 
complex orbit being given by the width of the corresponding peak in the Fourier transform, 
its exact value is hidden by the broadening of this peak. For the same reason, it is impossible 
to distinguish peaks closer one to each other than the "Fourier limit" as it is the case for 
the two orbits created at a bifurcation point. 

In this paper, we show how an excellent agreement between the quantum and the semi- 
classical properties in the vicinity of a bifurcation can be obtained using the harmonic 
inversion, a newly developed method, which is able to bypass the Fourier limitation |[T9|-pT|. 
The paper is divided as follows : in section |I| the essential properties of the hydrogen 
atom in a magnetic field are given, then in section |T| we present the usual semi-classical 
approach of the saddle-node bifurcation. In section IV we briefiy explain the harmonic 
inversion method, which is applied in section [V| to follow, from our quantum calculation, 
the properties of classical orbits through a saddle-node bifurcation. 



II. HYDROGEN ATOM IN MAGNETIC FIELD 

In atomic units, the Hamiltonian of an hydrogen atom in a magnetic field is given by 
(using cylindrical coordinates) : 

F = lp2 _ 1 + l^L, + ^7^' (1) 
2^ r 2' 8^ w 
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where 7 = B/Bq, with Bq = 2.35 x lO^T. Due to the rotational invariance around the 
^— axis of the Hamiltonian, is a good quantum number and we shall take = in what 
follows. 

The classical counterpart of this Hamiltonian has a scaling property, that is, if we define 
new variables by : 

< p = 7-V3p (2) 

we obtain a new Hamiltonian H given by : 

H = r^l^H = P- - i + ^. (3) 

' 2 f 8 ^ ' 

which does not depend anymore on 7. The classical dynamics of this Hamiltonian is entirely 
fixed by the scaled energy e given by : 

e = 7-2/3^. (4) 

All properties of the classical trajectories of the original Hamiltonian can be deduced from 
the scaled dynamics using the scaling transformation (^. For example, the action S of an 
orbit (i.e. / p.dq) is related to the reduced action S{e) (all quantities with "refer to the 
scaled Hamiltonian) by : 

S{E,^)=r'''S{e). (5) 

From the quantum point of view, this scaling introduces an effective h value, which is 
easily seen on the scaled Schrodinger equation. Hip = eip, for fixed scaled energy e : 

Thus, the effective h is given by 7^/^ and so at fixed value of the scaled energy e, the 
semi-classical limit is obtained when 7 tends to 0. 



III. SEMI-CLASSICAL APPROACH AND SADDLE-NODE BIFURCATION 



The semi-classical approximation of the photo-ionization cross-section (7(7, E) is obtained 
by a generalization of the Gutzwiller trace formula [II^P^I and has the following expression : 



cr(7, E) = crcoui(^) + o-osc(7, E). 



(7) 



o'couiiE) is a smooth background term and corresponds to the Coulombic cross-section that 
would be obtained in the absence of a magnetic field and, for energies close to the ionization 
threshold E = 0, aconi{E) is almost independent from E, so that in first approximation 
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o'Coni{E) ^ o"coui(0). The second term is the oscillatory part of the cross-section and is given 
as a sum other all orbits closing at the nucleus (and their repetitions) [|15,1C,2^ : 

aosc(7, E) = ^J2 Ml, E) sin (0^(7, E)) , (8) 
k 

uj being the photon frequency. Actually, looking at cross-section for energies close the 
ionization threshold (i.e. = 0), = {E ~ Ei)/h is almost given by the initial state energy 
—Ei and does not depend on the energy E. In the preceding formula, each amplitude 
Ak involves different classical quantities |2^, among which the matrix element mu of the 



monodromy matrix (i.e. the stability matrix restricted to deviations perpendicular to a 
closed orbit in the energy shell). More precisely, Ak has the following expression : 



Ak = 2{2nf/ysmeiSmefyU^^)ym{0f)^=, (9) 

\mi2\ 



where 6ij are respectively the initial and final angles (at the nucleus) of the trajectory with 
the axis z of the magnetic field. The function 3^m(6') depends only on the structure of the 
initial state and the polarization of the exciting light. Physically, it represents the angular 
distribution of the excited electron leaving the nucleus. Explicit expressions for excitation 



with TT-polarized light from the n = 1 (ground state) or n = 2 are the following [p^p3| , p5 



= |lsO) 


3^o(^) = 


-7r"^/22V^cos^ 


\^P,) = \2s0) 


3^o(^) = 


-{2n)-'/^2^e-^cos9 


\^,) = \2p0) 


3^o(^) = 


(27r)-i/22^e-^(4cos2^-l 


The phase 0^ is given by : 







(10) 



(Pk = Sk{-f,E)~^fik + ^, (11) 

where Sk is the action of the closed orbit. The index /i^ takes into account the different 
phase shifts occurring at conjugate points, at the nucleus and at the crossings of the 2;-axis. 
This index will be called Maslov index, even if it is not the usual one. More precisely, it is 
given by the following formula (for m = 0) ||15|| : 

/i = Z/Q + Z/l + Z/2 + Z/3, (12) 

where i^o is the number of conjugate points (mi2 = 0) along the closed orbit, ui is the number 
of points at which the velocity vanishes (only for self retracing orbits), z/2 is the number of 
crossings of the 2;-axis and z/3 is the number of times that an orbit leaves and reaches the 
nucleus. The case of the orbit along the field (i.e. the 2;-axis) is special and requires a slight 
modification (for details, see PJT5[]). 



The singularity in the classical equations of motion due to the divergence of the Coulomb 
potential at r = is regularized using the semi-parabolic coordinates [u = y/r + z, v = 
y/r — z), giving rise to the following effective classical Hamiltonian |p5| , p6| : 

h = lpl + \vl - + v^) + \v?v\u^ + v") = 2. (13) 



4 



One can show that, due to the potential shape, there is a deep analogy of the classical 
dynamics with the 4-disk problem [^-2^ (i.e. the dynamics of a point particle moving 



freely in the {u, v) plane and bouncing elastically off the walls made by four disks at the 
corners of a square). In the magnetized hydrogen atom, the role of the 4 disks are played 
by the 4 potential hills between the valleys along the u and v axis. Thus, the coding scheme 
of the 4-disk can be directly used : the (m, v) plane is divided in four parts delimited by the 
coordinates axis and each part is labeled with a number e {1, 2, 3, 4} (see Fig. ^). Then for 
each orbit, a sequence Si . . . s„ is built, corresponding to the sequence of bounces ||30| made 
by this orbit. At large positive scaled energy (above e = 0.33), it seems that this coding 
scheme works perfectly well . At lower scaled energy, some sequences of symbols become 
forbidden (pruning takes place); in the region studied in this paper, the coding scheme is 
still quite efficient. It allows us to derive few properties like the Maslov indices fik in a very 
simple way l^-H. 

Using scaling properties of the classical dynamics, the oscillating part of the cross-section 
can be rewritten as : 

t^osc(7,e) = ^Il7^^^f^fc(e) 

^ ^ (14) 

= ^ E j^/'Me) sin [r'''S,{e) - |/..(e) + ^) . 

The scaled amplitudes are functions of the scaled energy e only (except for the orbit 
along the field, which is varying like 7^"^^). The exact expression is directly derived from 
equation : 

ifc(e) = 2{27TfysmeiSmefyUO^)ym{0f) , ^ (15) 

|"^i2(e)| 



Looking at formula ([T^ ) one can see that, by fixing the scaled energy e and doing the 
Fourier transform of 7~^/^(j(7,e) with respect to the variable 7"^''^, one must obtain peaks 
at the scaled actions of the closed orbits : this is the well-known scaled spectroscopy [0]. 

However, the preceding approach fails in case of a saddle-node bifurcation at a given 
value ec of the scaled energy. At the bifurcation, two orbits are created, but with fhi2 
coefficients equal to and thus the scaled amplitude of the two orbits are infinite. This 
divergence can be regularized by a detailed study of the modifications induced in the trace 
formula by the bifurcation |lT8| , p5| . More precisely, the contribution to the oscillating part 
of the cross-section due to the orbits involved in the bifurcation becomes : 
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3ay/6 

2 ) 

X M ((f )^/^-/^(e. - e)) sin (7-^/^(^(e.)) - ^/.°) , (16) 

where Ai{z) is the Airy function [^, a and M are defined by local expansion of S and mi2 
near the bifurcation, that is : 



'-1 

\3/2 




S'(ec) ±(T (e - Cc 
±M (e - e,) 



1/2 (17) 
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where ± refers to the pair of orbits born above the bifurcation. At fixed value of the scaled 
energy e > ec and for very small values of 7 which corresponds to the semi-classical limit, 
using the asymptotic behaviour of Airy function |31|], one obtains : 



= ^7^/^io(e) sin (l-'/HS{e)_) - - 
+ ^I'^'Me) sin (l-'^HSieU) - |(/.° + 



(18) 



with Ao(e) is given by : 



io(e) = 2i27rfysme,sm9fymmyU0f) , . = (19) 

|M||e-eJi/2 



recovering thus the usual semi-classical contribution of the two closed orbits created at the 
bifurcation. 



Below the bifurcation and again using the asymptotic behaviour of Airy function |^ 
one obtains : 

= ^7^/^io(e) sin {j-'/'Sie,) - e--^^^-(--)^^^ (20) 

The contribution does not vanish, even if it is exponentially decreasing when going away 
from the bifurcation or going to the semi-classical limit 7 0. It has exactly the same 
functional dependence than for a usual closed orbit, if one allows the action to become 
complex S{e) = S{ec) — (ec — e)^^^ , which is nothing but the continuation of Eq. (|T^ 
across the bifurcation. In fact, it can be interpreted as the contribution of a "ghost" closed 
orbit, hving in a complexified phase space [ITSj^ . 



IV. HARMONIC INVERSION 

This method is a powerful tool to extract Fourier components (phase and amplitude) of 
a signal, with a better accuracy than one obtained by the usual Fourier transform, which 
is limited by the total length of the signal. Especially, this method is very well adapted to 
distinguish peaks closer than the Fourier resolution, provided that they are well separated 
from the other ones. A complete description of this method and some of its applications can 
be found in references Il2llj3^j33|| . 

Given a time signal c{t), known in the interval [0,T] and for which we assume the 
following expression : 



AT 

c(t) = ^ a^e-'--', (21) 

n=l 

where and Un are the unknown amplitudes and frequencies to be extracted from the 
signal. 
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Of course, these can be obtained using standard Fourier transform, but then the reso- 
lution is hmited by the interval length T, adding an artificial width equal to 2'k/T to each 
Fourier peak. Especially, for frequencies closer one to each other than 4:71 /T, their values 
does not correspond to the respective maxima in the Fourier spectrum, as one can see from 
the following simple example : 

where a and b are real numbers. The modulus square of the time limited Fourier transform 
(i.e. fiu) = /o^rftc(t)e'"*), is given by : 

|/(cu)p = a^T^ sinc^ {uj - ^+)| + b^T^ sinc^ (cu - cu_)^ 

T T T 

+ 2abT^ cos (0;+ — ^-)"^ sine {u — ^+)^ sine {u — ^-)^ (23) 

where the function sinc(x) is the usual sin(a;)/x. As expected, it is only a function of 

= 00 — ujq : 

\f{Co)\^ = a^T^ sinc^ {Co - ^5cu)| + b^T^ sinc^ {u + ^5cu)| 

T ^ 1 T ^ 1 T 
+ 2a6T^ cos Su— sine (tu — — sine (tu + -6uj) — (24) 

For T large enough (i.e. much larger than Tq = 2ti /5uj), the two peaks of the Fourier 
transform are well separated and we are able to extract from this Fourier transform the 
right values of the two frequencies uj^ = 5uj/2 and uj_ = —5uj/2. On the contrary, as one 
can see in Fig. 0, when T is of the order of 2ti/5uj, the effective positions of the peaks are 
less and less well defined. Even for T = Att/Suj (i.e. theoretical Fourier resolution two times 
as small as 5a;), the distance between the two peaks is enlarged by 20 percent. Furthermore, 
how this shift is distributed between the two peaks depends strongly of the relative amplitude 
a/b. For example if 6 is 5 times as small as a, then b is shifted by 85 percent of the total 
shift, whereas a is only shifted by 15 percent. 

It becomes worse and worse as the number of frequencies in the signal increases, which 
shows that the usual Fourier transform is of little interest to extract accurate information 
from a limited signal. On the contrary, the harmonic inversion can bypass this Fourier 
limitation. 

The main idea of the harmonic inversion is to construct an abstract evolution operator 
U(t) and an abstract initial state |$o) such that c{t) is given by ($o|^(^)|'^'o)5 where ( . ) 
is a complex symmetric inner product (i.e. {'ip\ is just the transpose of I?/'), without complex 
conjugation). The eigenvalues of this operator are the e"**^"*, showing that U is not an uni- 
tary operator {uj^ being generally complex). Then by diagonalizing its matrix representation 
in a suitable basis, it is possible to extract with a much higher accuracy than the Fourier 
limit these eigenvalues; the associated eigenvectors are then used to find the amplitudes a„. 
One possible basis is a Krylov basis : \n) = U{n6t)\^o), where St is a short time interval 
and n is an integer {0 <n < N, N defining the basis size). Thus, the matrix representation 



7 



is simply given by Ann' = {n\U{St)\n') = c(^{n + n' + l)St^. The basis vectors being not 
orthogonal {Bnn' = {n\n') = c({n + n')St]), one has to solve the symmetric generalized 



eigenvalues problem : 

A\(j)) = e-*"^*5|0) (25) 

However, this simplest basis choice is not satisfactory for efficient numerical purpose, 
because all coefficients in the matrices have the same order of magnitude, such that a 
eigenvector will have significant overlap with all basis states. Furthermore, all frequencies 
ujn are equally represented, which means that N has to be of the order of the total number 
of frequencies, which can be very large. For this reason, another basis, using the filtering 
properties of the Fourier transform is introduced. The choice of the basis is done by selecting 
an interval [cUinf , t^sup] in which one wants to extract frequencies. For frequencies Uj in this 
interval, one builds a new basis \ipj) = X^n From the signal c(t), only frequencies 
in the interval will give significant contributions to the basis vectors. Thus, one can take a 
much smaller basis, of the order of the number of frequencies in this interval. Moreover, the 
matrix structure of the evolution operator in this basis in essentially banded, that is matrix 
coefficients are decreasing rapidly when going away from the diagonal. The expressions 
of the coefficients and other details about the numerical resolution of the new generalized 



eigenvalues problem can be found in pl|,|32,33 



A first test of this method is to apply it to the preceding example. The results are shown 
in Fig. ^ For both graphs, the continuous line is modulus of the usual Fourier transform 
divided by T and the vertical segments are the amplitudes given by the harmonic inversion. 
Sizes of the Krylov basis from 2 to 50 have been used, although 2 would have been enough, 
to emphasize that even when using an over-complete basis (which would be the standard 
case, as the number of frequencies is unknown) the results given by the harmonic inversion 
are not affected. The upper graph has been made for a signal length equal to Tq = 2Tr/Suj, 
i.e. the theoretical Fourier resolution is the distance between the two peaks. The bottom 
graph is made for a signal length Tq/IO. From both pictures, it appears clearly that the 
harmonic inversion is much more efficient than the usual Fourier transform : for T = Tq, 
the Fourier transform still shows two peaks, but shifted by more than 20 percent from the 
exact positions and for T = Tq/IO, the two peaks have collapsed in a single peak, whose 
width is much larger than 6uj. 

On the contrary, values extracted by the harmonic inversion are very well localized, 
at the right frequencies. More precisely, the accuracy on different parameters (position and 
amplitude) are shown in table |. Of course, this example is rather simple, but even for signal 
length such that the Fourier resolution would be ten times as large as the peak separation 
and for relative amplitudes differing by two orders of magnitude, the accuracies are better 
than 10~^ on the frequencies and 10~^ on the amplitudes. For more complicated situations, 
the accuracy is expected to be worse than in the present example, but it will be still much 
better than what we would be able to extract from the Fourier transform. 
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V. RESULTS 



How one can apply the harmonic inversion in the case of the hydrogen atom in a magnetic 

corresponds to time t, the 



field is obvious if one compares equations ( pi]) and ([T^) : 7"^/^ 
scaled action Sk{t) to (minus) oOn and exp [—i'^jj.k + «f ) to a„. For negative scaled energy 
e, the function a (7, e) is a just a sum of Dirac delta functions : 



f^(7,e) 



7„-^/^) 



(26) 



n=0 



where 7^^^^ are the eigenvalues of the scaled quantum Hamiltonian (i.e. e being fixed, 
energy E quantization, because of equation (^, is equivalent to 7"^/^ quantization). /„ are 
the associated (squared) excitation matrix elements. 

That we explicitly know the functional dependency of c{t = 7"^^'^) can be used to simplify 
the implementation of the harmonic inversion. Indeed, the latter was written in the case of 
a general signal obtained either numerically or experimentally and known only at equally 
spaced values of time. On the contrary, in the present case, we are able to compute c{t) for 
any time t. A straightforward method would be to compute c{t) only on a grid. But because 
of 6 functions we would have to add an artificial width to each peak for numerical purpose, as 
it is done in ref [Q. The other way round is to modify a little bit the generalized eigenvalues 
problem ([25|), in which the matrix representation of the evolution operator U{t) was used, 
to obtain directly a generalized eigen-equation for the operator Q, defined by U{t) = e~*^*, 
which is very close to what was done in the original articles of Neuhauser |TP|J2^. If one 
writes the eigensystem in the following form : 



^1 

then the coefficients of A are given by : 

^(01,02) = 
(01 7^ 02) 



-iujBl 



(27) 



1 



n — HJ2 

; T 



T/2 /■T/2 

Jo 



2) 



^^¥^1-^^) r dtc{ty'^^ + (Pie'^^'t'^-'t"^ r dtc{ty'^' 

Jt/2 Jt/2 

rT/2 /■T/2 /■T 

dtc{t)e''^' dtc{t)te'"^' - dtc{t)e'"^^ 
/o Jo Jt/2 

+ iT(t)i r c{t)e''^' - Kpi r dtc{t)te''^' 
Jt/2 Jt/2 



(28) 



/T/2 JT/2 

where 0i and 02 are two frequencies taken in the interval [u;inf, u!sup] in which one wants to 
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extract the uon- The matrix elements of B have similar expressions : 

^(01,02) = 



7^ 



1 T V2) 



T/2 i'T/2 

Jo 
T, T/ ^ 1 (29) 

T/2 Jt/2 
T T T 



>) / dtc(t)e'"^' 

Jo 



' ' '2 '2 

Insertion of signal (^) in the preceding expressions is straightforward and leads us to replace 
all integrals with sums on the available eigenenergies Jn^^^- 

Using the preceding implementation of the harmonic inversion, we will focus on the 
so-called XI bifurcation which is of the saddle- node type : at the bifurcation energy 
ec = —0.11544216, two closed orbits are created, one stable and one instable. The important 
point is that these two orbits have classical actions very well separated from all other closed 
orbits, so that their evolution with the scaled energy can be easily followed (both on the 
experimental and theoretical points of view). 

All quantum properties (energy levels and excitation matrix elements) have been ob- 
tained by numerical diagonalization of sparse matrices, which are representations of the full 



Hamiltonian in sturmian basis (for further details see ref. ||3^). Thus, for each value of 
the scaled energy e, we have ~ 10000 eigenvalues, giving a length signal equal to ~ 120. e 
values are ranging from —0.15 to —0.07 with 0.001 step. For each scaled energy below ec, 
the complex scaled action and the amplitude of the ghost peak have been extracted, whereas 
for scaled energy above ec, actions of both real orbits created at the bifurcation and their 
amplitudes have been extracted. 

To emphasize the agreement with the semi-classical formula (|I^) (^), we have also com- 
puted, for each value of the scaled energy, the classical quantities (scaled action, Lyapunov 
exponent and Maslov indices). For scaled energy below ec, the ghost orbits are found by 
extending the classical equations of motion in complex plane, which also makes time to be 
complex. More precisely, for real closed orbits, one has to determine two parameters (i.e. 
the initial direction of the orbit and its time length), whereas for complex closed orbits 
one has to determine four parameters (i.e. the complex initial direction and the complex 
final time). However, integration for complex time requires to integrate along paths in the 
complex plane; it is far from obvious that the (complex) classical quantities are functions 
of the final time only. For simple systems, like a one- dimensional double well potential, 
one can show that the classical quantities are meromorphic functions of the final time, but 
to our knowledge, there is no exact proof in the case of the hydrogen atom in a magnetic 
field. Still, the potential being a polynomial function of the coordinates (m, v), the structure 
of the equations of motion is such that the classical quantities are locally analytic, so that 
small deviations from a given path lead to the same results. Thus, even if a complete study 
of all singularities and branch points in the complex plane would be needed to justify this 
approach, one restricts integration of classical equations of motion on straight lines in the 
complex time plane starting from the origin, the integration time being simply the length 
along the line. In that case, the usual Runge-Kutta method can still be used. For example, 
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the ghost orbit associated with the XI bifurcation is shown in Fig. H(top). The scaled en- 
ergy is e = —0.14. The continuous hne is the real part, that is in the plane (^Re(M), Re(u)^ , 

whereas the dotted line is the imaginary part (i.e. in the plane (lm{u), Im(f )^. On the bot- 
tom are shown the two real orbits (at scaled energy e = —0.11) created at the bifurcation. 
The dashed line corresponds to the S- orbit and the continuous line to the S+ orbit. The fact 
that the only difference between the two trajectories is the additional bounce (i.e. another 
conjugate point) in the (m > 0, i; < 0) part of the space for the longest orbit (5*+) explains 
that the Maslov indices associated with these orbits are /i_ = yUo and /i+ = /io + 1- If we 
use the coding scheme of the 4-disk problem discussed above, the two orbits involved in the 
bifurcation have respectively codes 124 and 1214 (see Fig. from which (using Eq. ([T^)) 
we calculate the Maslov indices /io = 8 and /io + 1 = 9. 

The first comparison is shown by Fig. | where scaled actions of the orbits are shown as 
functions of the scaled energy. Below the bifurcation (on the left), is shown the real part 
of the action and above the bifurcation (on the right), are shown the actions of the two 
orbits born at the bifurcation. The continuous lines are the classical calculation, whereas 
the circle are the values extracted from the quantum data. The agreement is excellent, even 
very close to the bifurcation point. This agreement is emphasized in Fig. |^, where on the 
left side, is shown the imaginary part of action of the ghost orbit, and, on the right side, is 
shown the difference AS" = (^2 — Si)/27r. Here again, the agreement is good, especially if 
one compares with the Fourier resolution which is ~ 0.0085. It shows clearly that it would 
have been impossible to extract the right behaviour of neither the width of the ghost peak 
nor the separation in the actions of the real orbits where it is the most interesting, that is 
as close as possible of the bifurcation. This clearly appears in Fig § where the usual Fourier 
transform, that is 



\F{e,S)\ 



(30) 



is drawn in the (e, S'/27r)-plane. It is obvious that one can not distinguish the two orbits for 
scaled energy e closer to the bifurcation than e = —0.085. Furthermore, even if for scaled 
energy greater than —0.085, the two real orbits seem to be far enough one from each other to 
be seen on the Fourier transform, we know from the simple example shown in section |V| that 
even for theoretical Fourier resolution twice better than the peaks separation, the extracted 
positions are still shifted by ~ 20 percent. This tells us that quantitative comparison with 
the classical dynamics would have been possible only for scaled energy larger than —0.07. 
But at that point, other peaks corresponding to a family of orbits start to overlap with these 
two peaks, so that it would be hardly impossible to extract any quantitative information. 
Also, it would be difficult, at present time, to have a Fourier resolution twice better only by 
increasing the signal length : as the density of states is proportional to 7"^^^, the number 
of states below a given 7"^/^ is proportional to 7"^''^, and thus to obtain a signal length 
twice as long, one would have to compute more than 40000 levels, which is at the edge of 
the present numerical possibilities (i.e. it would require 32 times as much CPU time and 8 
times as much memory). 

After having compared the frequencies, one can compare the amplitudes with the semi- 
classical predictions. The results are shown in Fig. ^ The continuous lines are the individual 
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semi-classical amplitudes as obtained from equation ([T5|), with the function correspond- 
ing to transition from the |2s0) initial state, with vr-polarized light (see equation ([To|) ). On 
both sides of the bifurcation point ec, two different regimes clearly appear. For |e — ed > 0.02, 
the agreement with the semi-classical prediction is quite good, even if not as good as for the 
actions. On the contrary for energies too close to the bifurcation point (i.e. |e — < 0.02), 
the discrepancy is large, showing thus that the semi-classical approximation is no more valid. 
This can be easily understood from equation (16), which shows that the semi-classical limit 
(i.e. 7~^/^ — >• -|-cxo) is in competition with the limit e — > ec- In other words, for fixed value 
of 7"^^^^, one goes out of the semi-classical limit when approaching the bifurcation. The 
break point corresponds to the argument of the Airy function of the order of one, which for 
7"^/'^ = 120 (the largest available value) gives e — ec ~ 0.01, in agreement with the results. 
For e > ec, the argument of the Airy function being directly related to the separation in ac- 
tion AS" (see equation (1^)), the preceding argument corresponds exactly to the fact that, for 
fixed scaled energy value, the phase shift 27r7~^/^AS^ has to be large enough (> 27r), so that 
the two orbits become distinguishable from the quantum point of view. A very simple test 
of the fact that we cannot use any asymptotic expansion of the Airy function, is to consider 
the amplitude exactly at the bifurcation point. In that case, the Airy function is only taken 
at the value 0, for any value of 7"^/^. More precisely, the term i^^^sin (^^~^^'^S{ec) — |yU° 
in the oscillating part of the cross-section (see equation (16)) becomes : 

= ^7^/' 2{27^r/y sin 9, sin Of yUe.)yUOf) (^^ 



2 



X \M\-'/^Az{0) sin [i-'/'S{ec) - (31) 

From this equation, one can see that it is possible to extract the amplitude from the 
quantum data using harmonic inversion, provided that the signal is multiplied by 7"^/^. 
For this purpose, 7000 levels have been computed at scaled energy e = ec. The amplitude 
extracted from this data is 2.94, in a very good agreement with the theoretical value 2.951 
given by equation (31). 

Having compared the modulus of the amplitude, we now turn to its phase, i.e. Maslov 
indices. As the Maslov indices occur only in the phases with a multiplicative factor tt/2, they 
can only be measured modulo 4. Fig. ^ shows the nice agreement between the Maslov indices 
extracted from harmonic inversion of the quantum signal and the semi-classical predictions. 
Indeed, from equation (|TBp, the Maslov indices of the two real orbits (e > ec) are respectively 
fj.0 and /io + 1 (with /iq = 8, i.e. modulo 4) whereas it is /iq for the ghost orbit (e < ec) 
which is exactly what is depicted by the figure. Exactly at the bifurcation point, with the 
modified harmonic inversion described in the previous paragraph, we obtain a Maslov index 
of 0.05, again in perfect agreement with the theoretical value /io = (modulo 4). 

Even if the agreement between the different quantities is excellent, it is not as exceptional 



as in the simple example (|2^) . The possible explanation comes from the fact that one 
assumes the expression for the signal to be a sum of exponentials with constant coefficients. 
This is clearly broken in the case of the semi-classical approximation because of the smooth 
contribution in either the density of state or the oscillator strength (i.e. Thomas-Fermi 
like terms) and also because of the remaining terms in the asymptotic expansion in h. All 
this will contribute to (slowly) varying coefficients, which will make the harmonic inversion 
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method less efficient. This is clearly emphasized by the behaviour of the amplitudes for 
scaled energy too close to the bifurcation for which one cannot use anymore the asymptotic 
expansion of the Airy function. 

VI. CONCLUSIONS 

In summary, we have shown that, using the harmonic inversion method, the properties 
(actions, Maslov indices and stabilities) of the classical orbits involved in the photo-ionization 
cross-section of the hydrogen atom in a magnetic field can be extracted with a much better 
accuracy than the usual Fourier transform, and this even for scaled energies close to a bifur- 
cation point. Below the bifurcation, the contribution of a ghost orbit has been emphasized 
by showing that the behaviour of the imaginary part of its action is in perfect agreement 
with the classical predictions. For scaled energy above the bifurcation, we have been able 
to distinguish the contributions of the two orbits created at the bifurcation and we have 
also shown the perfect agreement with the semi-classical prediction for the frequencies, am- 
plitudes and phases of the modulations. We have also emphasized the non-semiclassical 
behaviour of the amplitudes too close to the bifurcation point. 

Laboratoire Kastler Brossel is laboratoire de I'Universite Pierre ct Marie Curie et de 
I'Ecole Normale Superieure, unite mixte de recherche 8552 du CNRS. CPU time on a Cray 
C90 computer has been provided by IDRIS. 
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FIGURES 



1.0 




-1.0 ' ' ' ' ' ' ' ' ' ' ' 

1 2 3 4 5 



FIG. 1. Effective frequencies extracted from the signal c{t) = ex.p{i6ojt/2) + exp{—iSujt/2) 
using finite Fourier transform, i.e. dt c{t) exp iCot, as functions of T. For large values of T (much 
larger than Tq = 2tt/5u)) we recover the right values of cD, that is ^6lo/2 (dashed lines). On the 
contrary, for T of the order of Tq, the frequencies given by the Fourier transform are substantially 
shifted from the correct values. For example, for T = 2Tq, the distance between the two peaks is 
enlarged by 20 percent. 
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FIG. 2. Comparison between the usual Fourier transform and the harmonic inversion technique 
for the signal c(t) = exp (i5u;t/2) +exp {—i5ujt/2). For both graphs, the continuous line is the usual 
Fourier transform \F[T)\/T (T is the length of the available signal c(t)), whereas the vertical lines 
arc at the frequencies given by the harmonic inversion, their heights corresponding to the associated 
amplitude (also obtained by harmonic inversion). The upper graph is made for T = Tq = In/dco, 
which means that the Fourier resolution is of the order of the peak separation. Again, we recover 
the fact that the effective positions given by the Fourier transform are quite shifted from the 
exact values and the amplitudes are also not correct. On the contrary, all values (amplitudes and 
positions) given by the harmonic inversion are in a perfect agreement. The lower graph is made 
for T = To/ 10, i.e. the Fourier resolution is 10 times as large as the peak separation, so that the 
Fourier transform is almost constant in the interval [—1, 1] and thus gives no information on the 
position of the peaks. In contrast, the harmonic inversion method is still working, with a very good 
accuracy (better than 10~^). 
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FIG. 3. Closed trajectories of the electron of an hydrogen atom in a magnetic field in the {u, v) 
(semi-parabolic coordinates) configuration space. On the top is shown the complex ghost orbit 
associated with the XI peak experimentally observed in the scaled spectroscopy of the hydrogen 
atom in a magnetic field at scaled energy e = —0.14 [^. The continuous line is the real part of 
the trajectory, whereas the dotted line is the imaginary part. Eventually, this complex trajectory 
collapses with its complex conjugate in a saddle-node bifurcation at €c = —0.11544216, from which 
two real orbits are created, shown on the bottom at scaled energy e = —0.11. The dashed line 
corresponds to the S- orbit created at the bifurcation, whereas the continuous line corresponds to 
the S+ orbit. From this plot, one can easily deduce the code of each orbit, that is the sequence 
of bounces in the different quarters of space (1,2,3,4) defined by the coordinate axis : for the 
code is 124 and 1214 for (one additional bounce on the (u > 0,f < 0) part of the plane). 
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FIG. 4. Real part of the scaled actions as function of the scaled energy e. The continuous 
lines are the classical results, whereas circles are extracted from the quantum dynamics using 
harmonic inversion. One can follow the two real orbits born at the saddle-node bifurcation, even 
for separation much smaller than the standard Fourier limitation 0.01). 
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FIG. 5. Usual scaled Fourier transform (see equation (|30|)) of the photo-excitation cross-section 
of the hydrogen atom in a magnetic field, in the (e, S'/27r)-plane. Below the bifurcation 
€c = —0.11544216, one clearly sees the contribution of the ghost orbit. However, it would be 
very hard to extract the relevant information about the complex orbit, namely the imaginary part 
of the action, from the width of the peak. Especially, it should vanish like (ec — e)'^/^ near the 
bifurcation, which is not the case in the figure. Actually, it is dominated by the broadening due 
to the limited range of the Fourier transform. Above the bifurcation, the actions of the two real 
orbits created at the bifurcation are separated by the same power law. But the separation cannot 
be observed on this plot below e ~ —0.085, which is already too far from the bifurcation point to 
observe this power law. 
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FIG. 6. On the left is the imaginary part of the scaled action of the ghost orbit (below the 
bifurcation) and on the right is the separation in action between the two real orbits above the 
bifurcation. For both graphs, we include levels up to jmlL'^ = 120, so that the theoretical Fourier 
limitation would be 8.5 x 10~^, but (see text), the actual resolution is at least twice as large. 
On the contrary, as it appears clearly on both figures, the limitations of the harmonic inversion 
are much smaller, which allows us to emphasize the good agreement between the semi-classical 
predictions and the exact quantum results. 
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FIG. 7. Comparison of quantum amplitudes ^fc(e) (see equations. ( [1^ ) (^)) extracted using 
the harmonic inversion, with the semi-classical prediction. The continuous lines are the amplitudes 
calculated from the classical quantities (stability, initial and final angle). Diamond (ghost orbit), 
crosses and circles (real orbits) are the amplitudes extracted from the exact quantum results. 
On both sides of the bifurcation, two regimes appear. For scaled energy away enough from the 
bifurcation (|e — ed > 0.02) the semi-classical predictions and the quantum results agree well. On 
the contrary for scaled energy too close to the bifurcation point, the deviation is large. 
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FIG. 8. Maslov indices (see equation (|l^ ) in text) associated with each orbit extracted from 
the exact quantum results. Above the bifurcation point ec, we recover the semi-classical prediction 
that Maslov indices of the two real orbits created in a saddle-node bifurcation are /io and /io + 1, 
where /io is also the Maslov index of the ghost orbit. Furthermore, this value = (modulo 4) 
agrees with the one obtained from Eq. ((|l^)), giving fiQ = 8. The symbol (+) correspond to the 
yUo value extracted exactly at the bifurcation point, using a specific harmonic inversion technique 
(see text), also in perfect agreement with the theoretical predictions. 
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TABLES 



TABLE I. Accuracy on parameters extracted from c{t) = aoe-^^'^°+^'^/^^* + 6oe-*('^o-'^'^/2)* 
using harmonic inversion, for signal length T = IOTq, T = Tq and T = Tq/IO, where Tq = 217/600. 
This table emphasizes the excellent accuracy of the method. For T > Tq, it is of the order of the 
machine precision. For signal length ten times as short as the theoretical Fourier resolution, the 
accuracy is still of the order of 10"^ on the frequency positions and 10"^ on the amplitudes. Even 
for very peculiar situations (one peak 100 times as small as the other one and signal length 10 
times as short as the theoretical Fourier resolution), the accuracy is still very good. 
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